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ABSTRACT 

This paper presents the heat transfer from a spherical balloon type 
shadow shield to a spherical cryogen storage tank under space conditions. 
The analysis covers these design factors; 

1, Combined skin eonductanci and internal radiation for the shield 
High emittance bands on the iniide and outside of the shield 

3. Optimum width of these bands 

4. Timpirature-varying-abiorptanci and imlttance of the sMeld’s 

surface 
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The priiint analyiii § istends that of Nicholi [8] who imposed constant 
emittance. And it generaliiii the work of Jones and Barry [2, 8] who as- 
sumed Isothermal radiating lones. We obtalnid the following rtiultii 
1, The high emittance bands and timperaturi varying tmittance com- 


bine to lower shield timperatures significantly, Banding allows more heat 
loss to space, Timpirtturt varying imittanci narrows the spread in tem- 
perature in going from the sum side to the shadow side of the shitldi and it 


irs the shiild timperaturi level, 

i. Joint inside and outside banding is much more effectivi than outside 
alone, Optimum inside and outside banding reduces boil- off heat 
flux by a factor of 4, Outside banding alone gives at most a 12% reduetlon. 

3, Optimum high emittance band width is a weak function of shield-tank 
separation. The optimum band width is 18® at a separation of one tank 
diameter, - (tank and shield ^ameter are the same in this analysis). 
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SUMMARY 

This report presents the heat transfer from a spherical ballon type 
shadow shield to a spherical cryogen storage tank. 

To find the heat transfer it is necessary to find the temperature pro- 
file around the shield. We do this by solving the heat equation for combined 
radiation and conduction heating of elements of the shield. The equation ac- 
counts for the presence of a high emittance band about the shield; and it in- 
cludes the effect of temperature-varying-emittance. Because of the form 
of the boundary condition, the equation required an iterative machine solu- 
tion of the Runge-Kutta type. For small conductance, the temperature is 
nearly constant on the shadow side of the shield. This constant tempera- 
ture occurs even though the emittance changes abruptly in crossing the high 
emittance band. The solution shows that the effect of the high emittance 
band is to lower the shield’s temperature level. Temperature-varying- 
emittance not only lowers shield temperature; it also narrows the spread 
in temperature in going from the shield’s sun side to its shadow side. 

Using the temperature solution and radiation view factors, found sep- 
arately, we predict the heat flux to the spherical storage tank. In this 
analysis the shield and tank diameters are equal. We present figures 
which show the heat flux as a function of width of the high emittance band 
and shield-tank separation. There is an optimum band width. This opti- 
mum band width is a weak function of separation distance, -varying from 
15° to 25° as the separation varies from 0 to 4 tank diameters. The heat 
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transfer calculation further shows that joint inside and outside banding is 
much more effective than outside banding alone. Two side banding reduces 
the heat flux by a factor of 4 whereas one side banding gives a reduction of 
only 12%, 

INTRODUCTION 

The aim of this paper is to facilitate the design of cryogenic storage 
systems for spacecraft. Specifically this report gives typical design 
curves for the heat transfer from a spherical shadow shield to a spherical 
fuel tank (of equal radius). Curves like these can be useful in the initial 
design stages in estimating fuel requirements, and in trade-off studies 
considering various shielding methods. 

Interest in long term storage stems from planned interplanetary mis- 
sions which call for cryogenic fuels such as liquid hydrogen (LHg) and liq- 
uid oxygen (LOX), For example, a Mars mission involves a round trip 
time of over 400 days. Over such a long period, poor thermal design of the 
cryogenic storage system could lead to excessive fuel boiloff. On the other 
hand optimum thermal design gives maximum payload and efficient lift off. 
During the planet transfer portion of the flight, the largest external 
heat source will be the solar heat fliix, Smolak, Knoll, and Wallner [1] 
have shown that shadow shields can greatly reudce solar heating under such 
conditions. Inflated shadow shields [2, 3] promise low structure weight 
and, therefore, high shielding efficiency, - in terms of boiloff weight to 
shield weight. 

However, there are other factors which enter into the choice of a ther- 
mal control system. Such factors may be the added fuel needed for vehicle 
orientation required with some systems, the presence of on-board heat 
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sources, or the presence of several external heat sources. For example, 
if solar and planetary heat fluxes are present, the designer may consider 
encapsulation of the cryogen module within a multi-layer foil. In any event, 
the system for a specific mission must be chosen from several known 
shielding methods. The results of this report simplify the consideration of 
spherical balloon shields when making such a choice. 

Specifically, this analysis finds the temperature profile around a thin- 
walled, hollow sphere subject to solar heating. And then, using this pro- 
file, the analysis predicts the heat transfer to a spherical storage tank in 
the shield's shadow. The temperatures result from a solution of the in- 
homogenous heat equation. The equation accounts for the distribution of the 
incident heat about the sphere by skin conduction and internal radiation. 
Furthermore, the equation allows for variable shield surface emittance. 

The shield emittance may vary with skin temperature, e = e(T); and it may 
vary in a stepwise manner, - viz. , banding. Knoll, et al. , [4], in their 
study of flat plate shields have shown that well placed high emittance bands 
can reduce shield temperature and thereby increase shield efficiency. 
Herein, we find the effect of high emittance zones (blackened bands) placed 
on the inside and the outside surface of the shield, on its shadow side. 

These bands might result from painting a high emittance coating on an 
otherwise highly reflective material like aluminized mylar. If both the in- 
side and the outside are painted, we call that simple banding. If just the 
outside is painted, we call that one-side banding. 

This paper extends the work of Jones and Barry [2, 3] who also con- 
sidered spherical, banded shadow shields; and it extends the work of 
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Nichols [5] who first found the temperature profile around a thin- walled 
hollow sphere subject to solar heating. Jones and Barry neglected conduc- 
tion along the skin and, instead, assumed isothermal radiating zones for 
use in an electrical analogue solution. Nichols included skin conduction 
but he imposed constant inside and outside surface properties. Thus, his 
analysis could not include banding and temperature variable emittance, 

LIST OF SYMBOLS 


A 

b 


Fij 

H 

I 


lin 


^sol 

h2 

k 

L 


Q 

R 

S 

T 


area 

thickness of shield’s skin 

view factor of surface j for flux leaving surface i 

total radiant flux leaving a surface 

total radiant flux arriving at a surface 

internal impinging flux 

solar impinging flux ® cos 9 

impinging flux due to presence of tank 

thermal conductivity 

vector distance between differential radiating areas 

heat flux 

radius of sphere 

shield-tank separation distance 

temperature 

sealing temperature, temperature to which the shield would come if 
k = 00 , = G. = el thus 4€^(TT^ = Q!*# 
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Z(0) unit step function, Z{9) =1, 0 :< 9 < tt/2; Z{9) ~ 0, tt/2 < 0 < it 
a absorptance 

Q!g solar absorptance 

Q!j sin 0 ’ d0’ 

i3 ratio of inside to outside emittance 

e emittance 

e emittance at 

rj stretched spatial coordinate 

11 dimensionless skin conductance parameter 

# parallel solar flux vector 

angle between surface normal and vector L 
dimensionless temperature 
cr Stefan- Boltzmann constant 

0 spatial coordinate - spherical colatitude 

Subscripts: 

1 refers to shield 

2 refers to tank 

i inside 

o outside 

TEMPERATURE DISTRIBUTION ANALYSIS 
s This section gives the main points of the derivation and solution of 

the heat equation. For more detail, see appendix A. 


average absorptance over inside area of shield a 


-r 

2 Jo 



6 


Figure 1 shows the physical arrangement we analyze. Parallel solar 
heat flux arrives along the line of centers of the shield and the tank. Fig- 
ure 1 also shows a control volume containing a typical portion of shield 
skin. The steady -state energy balance for the control volume is: 


(Net conduction flux) == < 


leaving radiationj 

arriving radiation) 


\ 

j 

outside + inside J 

outside + inside J 


The arrangement in figure 1 has axial symmetry; further we assume that 
there is no radial temperature drop through the skin of the shield. Thus, 
there is conduction only in the 0 direction. Hence, the energy balance 
leads to a heat equation in the form 


kb 


sineSu (H„ ^ Hj) - + Ij2> 


sin 0 ^ 


60 


( 1 ) 


For definitions, see the List of Symbols. In particular 1^2 to 

the shield due to the presence of the tank. In full, 


'12 








oT^ + (1 - a 


0 ) 112 ) 


dAj^dA2 


dA. 


^ ^dAgdAj ^2 


( 2 ) 


For the problem at hand, 1^2 ts negligible because: 

A 


1. For a cryogenic tank at Earth orbit, 0 T 2 « 

2. For radiation view factors between spheres 

^dAjdA2^dA2dAj^ ^sol 


at Earth's orbit. 



7 


This is fortunate for in the absence of 1^2 equation (1) becomes an ordi- 
nary differential equation for T instead of an integro-differential 
equation, - viz. 


kb 


d 


dT^ 


R sin 8 




(3) 


After inserting the definitions of the incoming and outgoing fluxes: 

Ho = cX (1 - “Aol’l 

Hj = ^ (1 . V 

Isoi = ^ 

equation (3) becomes 

kb _ 

de 


(4) 


sin 9 


^ fsin 9 

d0 


-^= (e^ + €j)aT^ - « J. 


- cos 9 
in s 


(5) 


The full derivation (appendix A) shows that for a simply banded sphere, 
(i. e. , or ^ = 1) 


4oiI. 


in 


a # 
s 


Here oi is the average inside surface absorptance 
sin 9 d0 


a 


1 ft 

= ± \ a. SI] 

2 Jo ^ 


( 6 ) 


(7) 


Notice that depends upon position not only because of banding but also 
because it depends on temperature which is a function of position. The 
final form of the heat equation is 


kb 


A /sin 0 


R^ sin 9 ^ 


69 


\ A /“-V^ 

-1= + e.)^T^ - — 1-^1- cos 0Z(0) (8) 

J \Q'/\4 / 
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We follow Nichols and use as a scaling temperature, T^= 

This is the temperature to which the shield would come if it had infinite 
conductance, and no high emittance band. For thermal variation of emit- 
tance and absorptance of unblackened areas of the shield, we assume: 

T and ^ ~ ^ 




the heat equation becomes 

-JL- A. fsin e — - 4 cos d Z(0) 
sin 0 d0 y d0 y \ 2/ 

The boundary conditions are: 

d£=o e = ®=o, e = i 

d0 d0 


(10a) 


(10b) 


Equations 10(a), and (b) require an iterative solution since the con- 
stant W can only be found from the solution. Because of the form of the 
boundary conditions (10b) we found a numerical solution using a two-point 
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Runge-Kutta technique. The solutions of equations 10(a), and (b), are 
shown in figure 2. With the temperature profile known it is possible to cal- 
culate the heat transfer to the tank. 

Heat Transfer 

In finding the heat transfer, we assume the tank to be isothermal 
(i.e. , Tg independent of 0g) '^2 ^ case, the net 

heat transfer from the shield (body one) to the tank (body two) is 


Qj2" 


fr 


^ 2 ^ 12 ^ 1 1 


(1 - <aj^)(l - ®2^^12^21 


'shield 

where the radiation view factor is given by: 


Fi2 “ 


/ cos 2 cos ^ 
irL^ 


^2 ^2 ^1 


( 12 ) 


Equation (12) gives the heat transfer as a function of shield emittance and 
shield tank separation. The emittance variation €2 = €-^{ 9 ) accounts for the 
high and low emitting portions of the shield, - due to banding and tempera- 
ture variation. The shield-tank separation appears in the view factor. We 
present the view factors separately in [6]. 

DISCUSSION OF RESULTS 


Figure 2 shows the effect of skin conductance, ju, on the dimensionless 
temperature distribution around the shield. Figure 2(a) is for constant un- 
blackened surface emittance. In finding these temperatures we used a rep- 
resentative value for emittance, such as that of aluminized mylar c = 0. 045 
over the unblackened portion. We used an emittance of 1 for the 15° 
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blackened band. Figure 2(b) shows the temperature profile for the case 
where emittance varies with temperature. Here we used e^= 0. 0585. The 
relationship between the cases of constant and variable emittance is shown 
in appendix B. 

Figures 2(a) and (b) both show that the temperature profile for ii = 0.1 

is only slightly different from the profile for ii = 0. For materials and 

-3 -5 

skin thicknesses of technical interest, ix is between 10 . and 10 . There- 

fore for the rest of this paper, we discuss only the temperatures and heat 
transfer for zero shield conductance. This has the advantage that for 
ju = 0, the heat equation becbmes algebraic arid its solution is much easier. 

Finally, note that both figures 2(a) and (b) show a uniform temperature 
on the shadow (tank), side of the shield, for p. = 0. Uniform temperature 
occurs even though the emittance changes abruptly at 9 = 7r/2 and 

0 =- + 9,. This supports the Jones and Barry assumption [2, 3] of iso- 
2 ° 

thermal radiating zones. However, the present solution has the advantage 
of giving the band temperatures explicitly. 

Figures 3 and 4 compare the temperature profile for five analytical 
cases. As shows in the legends of the figures, these five cases are com- 
binations of one side or two side banding with and without temperature vari- 
able unbanded emittapce. The cases as defined in figure 3 and 4 will be 
used throughout this discussion. Case I, the standard, is the Nichols pro- 
file [5] for constant abSorptance and emittance (i.e,, neither banding nor 
e=e(T)). 
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In figure 3, profile II shows that banding and temperature variable 
emlttance lower shieJd dimensionless temperature, xj/, relative to case I. 

The decrease in dimensionless temperature, results directly from 
banding. In addition, the scaling temperature, T^, for case Il ls lower 
than for case I. This results from the increase in emittance with increas- 
ing temperature, Thus the dimensional temperature T = will be 
lower in case II than in case I. 

Further, profiles II and III, in figure 3, are continuous even though 
the emittance is discontinuous crossing the blackened band. For /x = 0, 
the temperatures could be discontinuous. Here they are not, for the fol- 
lowing reasons. Outside blackening increases the ability to lose heat to 
space; but this greater heat loss is just balanced by greater heat gain from 
internal blackening. Thus the temperatures are continuous. Profile IV, for 
outside black band only, strengthens this argument; because it shows a 
sharp drop in temperature at the blackened band. 

Figure 4 shows hov/ temperature varying emittance alters profile I 
without targetting. Notice that the spread in temperature from the sun side 
to the tank side is smaller when e - <£(T) . 

Figure 5 shows how the heat transfer from shield to tank varies with 
high emittance band width For clarity, we present these curves in 
dimensionless form. Shield temperature goes down as grows, but the 
emissive power of the shield increases with These two opposing ef- 
fects give rise to an optimum band width (minimum ^ particular 

S value). However, the view factor, and thereby the effective emittance, 
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depends upon shield tank separation. So the optimum band width also de- 
pends upon separation distance. The optimum band width varies from 12° 
for a separation of 8R to 21° for zero separation. At the optimum band 
width, targetting can reduce the heat transfer by as much as 75% (at zero 
separation) . This reduction decreases with increasing separation. The 
reduction in heat transfer is 60% at a separation S = 8R, 

Figure 5 further shows that, increases in bandwidth beyond optimum 
lead to heat transfer almost twice the unbanded value when 0^ = 5t/2. 

Figure 5 shows how the boil-off heat flux varies with high emittance 
band width and shield tank separation distance. These curves are repre- 
sentative of design curves which might be used in choosing between shield- 
ing methods. Since the curves require particular emittance values (e^, e^) 
they are not general. 

The ordinate of figure 5 is related to the solar heat flux as follows 



It is clear from figure 5 tliat one side targetting is far less effective than 
two side targetting in reducing heat flux. 

CONCLUSIONS 

1. Use of a high- emittance band on spherical shadow shields decreases 
heat flux to a spherical shielded body by as much as a factor of 4. For 
cases considered here, the optimum high emittance band width is between 
12° and 25° for separation of zero to eight tank radii. 
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2. Banding and temperature variable emittance both decrease shield 
temperature level, and hence reduce the heat flux. Further, variable emit- 
tance narrows the spread in temperature about the shield. 

3. Having high emittance bands inside and outside is much more effec- 
tive in reducing shield temperature than outside’ banding alone. 

4. For small conductance < 10”^), shield temperature is nearly 
constant on the shadow side. This supports the assumption used in the de- 
sign method of Jones and Barry [2, 3]. However, the present calculations 
are more simple and direct than those of references 2 and 3 . 

APPENDIX A 

DERIVATION OF HEAT EQUATION 


This section presents the full derivation of the heat equation, equa- 
tion (10) of the main text. 

From the energy balance 


kb 


^ (sine^" 


sin 9 ^ 


69 




(Al) 


The symbol H stands for the total amount of flux leaving a surface, (i.e., 
emitted plus reflected). For example 

Hi = ^ (1 - “i)IinS + d - ^ U ' “)li2 

The symbol 1^2 total flux impinge on the shield (body 1) because of 

presence of the tank (body 2) 


^12 


/I 


^ 2^'^2 ■ “ 2 ^ 


/ 

+ (1 - “i)Ii2]*'dAjdA2 ‘^^ij^'dAjClAj ‘'^2 
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The inside impinging flux is 

From Nichols [5] 

dAJ = 2ttV? sin 0’d9'; 

4ttR^ 

and 


(A3) 


lin(0) = - C [e^crT'^(0’) + (1 “ Q!i)lin^0’)]sin 0’ d9 ' 

2Jq 


(A4) 


Thus is independent of 0 , 


a = i / a. si 

2 J ^ 


sin 0’ d0’ + (1 - 


in 


(A5) 


sin 0’ d0' 


Return to equation (Al) and consider 1^2' For a cryogenic storage tank 


«L.i- So 


'^'^2 ^sol- spheres (c^crT Fd^^dA^^dA^dAg ^sol 

neglect equation (Al). Then (Al) becomes 


kb_ ^ 
d0 


[sin 0 


«2 . ^ 
R sin 0 


— V ('n + - “A 

d0 y ° ^ ■ ■ 


. - Q! I , 

in s sol 


(A6) 


Use (A6) to replace e^crT^ under the integral in (A5) 
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— sin 0 ’ — + Q!.I. 


+ Z(0)Q! $ cos 0’ sin 0’ d0’ (A 7) 

s 


Once we specify the variation of and over the shiold, the integral 
can be performed to make (A6) an ordinary equation. Suppose that € is a 
linear function of temperature; and let = |3 Cq. Then 


Gf + 1 + ^ 


Then 


i3ofg$ 


/3 = 1, 4Q!I.^ = Q!g® 


As a result (A6) becomes 


’ kb d 


sin 0 


sin 0 


(1 + ^)€QaT 


4 /“i 


- Q!„<& COS 0 

s s 


Using 


«3«; M 




sin 0 d0 


iL_ i- fsin 0 ^\= (1 + - (—¥ - ^ CO® ^ 
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for 


0 


0 < - 


and 


IJ- 


sin 9 de 


[sin 0 ® 
d0 


= (1 + 



for 


IT 


< 9 ^ 7T 


(AlOb) 


In the main text we assume /3 = 1 and 0!^= e^. The assumption i3 = 1 
implies that the outside surface is prepared exactly like the inside surface. 
For example, doubly aluminized mylar painted black inside and outside on 
the shadow side (viz. jr/2 9 < ir/2 + 9y^. K in addition to J3 ~ 1, we assume 
e = (e^/T^)T = it follows that 


fi d 
sin 0 d0 


^sin 0 


d0 


2i//^ - 



4' - 4(z(9)] 


COS 0 


(All) 


For the case of outside targetting alone, equation (A8) does not have 
such a simple form for the ratio (ej/(ej + e_)) does not remain constant. 
Suppose, for instance, that c. - = constant except over the targetted 

band where c,. = 1. Then 


a.I! 



! 

- 


2(1 + €j:) + sin 0|^(1 " €j) 

... kb 


^^i /dT\ 

4(1 + ej 

4r2 

\d0 / 

V2- (1+^1>\*A/2+ 

— -J 


^ cos0 


4R" 




2Ci 


('d'rX 


1 + Cj \d0 


/ [(7r/2)+0, - 


a 0 
s 

8 


(A12) 
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The symbol (dT/d 0 )^^ 2 - stands for the derivative approaching 7 t/ 2 from 
the left; and (c^T/d0)^y2+ i«ipli®s the value approaching from the right. 
However for k 5^ 0 but continuous, equation (A6) implies that [8] 



Even if the integrand has jumps, the right hand side must be continuous; 
consequently the left hand side must be continuous. Then (A 12 ) can be 
placed in (A6) to obtain 


11 d 
sin e d0 


X 


^2 + M 



.,.4 




(1 + e.) 


“00 


cos 9 


2(1 + €j) + (1 “ G.)sin 6 -^ 

- 4 Z(e) cos 9 (A 13 ) 


7 t/ 2 - 


"l^/[(7r/2)+0j 


This is the heat equation where only the outside is blackened. It requires 
an iterative solution. 


APPENDIX B 

I 

EVALUATION OF T^ FOR e A LINEAR FUNCTIONS OF TEMPERATURE 
We assume as an example, a variation of emittance with temperature 
of the form e - (^oo/Too)'^ for all temperature ranges of interest. This 
form was motivated by the form given by Siegel and Howell [6] for tempera- 
tures around room temperature] g = (g 2 I 73 / 273 )T. We used this latter 
equation to find for case H as follows: 
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For € 2 f 73 we used the handbook value for aluminum foil, e = 0. 045. This 
led to 

e^= 0.0585 
T^= 355‘^K 

For case in we used e = 0. 045 and 379° K. 

APPENDIX C 

MAGNITUDE OF HEAT FLUX TO THE SHIELD DUE TO THE TANK 
The analysis of the temperature profile neglected the heat flux to the 

shield due to the presence of the tank. The terms in equation (10) are all of 

4 o 

magnitude eoT^. Since the tank is at cryogenic temperatures (say 20 K), 



and direct emission from tank to shield is negligible. The tank reflects 
heat from the shield back to the shield. This heat is 

(^^l)^ 12^21^^ ~ 

^ " ^12^21^^ ” “2^^^ " ^1^ 
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Thus the relative size of this term is 

1 - Fi2*'2i(1 - “2X1 - "1) 

There are four cases to consider. These are the combination of the tank 

reflecting heat from the shiny or blackened portions of the shield back to the 

shiny or blackened portions of the shield. In each of the four combinations, 

-3 

the relative size of the flux was less than 10 for a shield-tank separation 
of one tank diameter. Thus there was small error in neglecting the pres- 
ence of the tank in calculating the shield temperature profile. 
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Figure 1. - Spherical shield analysis geometry. 
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Case 

I No blackening (Nichols [53) e = const. 

II Black band, inside and out 

plus tempeaature variable 
emittance, E 

Tir Black band, inside and out 
e independent of ijf 
IV Black band outside only, 

£ independent of if/ 



Figure 3. - Effect of targetting on dimensionless temperature 
profile (ji = 0. 9b= 15°). 


I Constant emittance (Nichols [53) 
V Emittance given by £ = e^iJj 

No blackening In either case 


4 



Figure 4. - Effect of temperature varying emittance (dt)) on 
dimensionless temperature profile (ji = 0). 
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